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The leading twist transverse momentum dependent parton distributions (TMDs) 

are studied in a light-cone description of the nucleon where the Fock expansion is 

truncated to consider only valence quarks. General analytic expressions are derived 

in terms of the six amplitudes needed to describe the three-quark sector of the nucleon 

light-cone wave function. Numerical calculations for the T-even TMDs are presented 

in a light-cone constituent quark model, and the role of the so-called pretzelosity is 

J> ■ investigated to produce a nonspherical shape of the nucleon. 
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In recent years much work has been devoted to study semi-inclusive deep inelastic scat 
tering (SIDIS) and Drell-Yan (DY) dilepton production as powerful tools to understand the 
nucleon structure. According to the factorization theorem that separates the coherent long 
distance interactions between projectile and target from the incoherent short distance inter- 
actions , the physical observables of such processes can be expressed as convolution 
of hard partonic scattering cross sections, parton distribution functions (PDFs) and parton 
fragmentation functions (FFs) . With respect to the usual inclusive deep inelastic 
scattering (DIS) where PDFs only depend on the longitudinal momentum fraction carried 
by the parton, now PDFs, as well as FFs, also depend on the transverse momentum. At 
leading twist there are eight transverse momentum dependent PDFs (TMDs) 

lid, 



three of them surviving when integrated over the transverse momentum and giving rise to 
the familiar parton density, helicity and transversity distributions. 

TMDs contain rich and direct three-dimensional information about the internal dynamics 
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of the nucleoli. In particular, they can help in understanding how the nucleon spin origi- 
nates from the quark spins and the orbital angular momentum of quarks and gluons. How- 
ever, being typical nonperturbative quantities, TMDs are not directly calculable in quantum 
chromodynamics (QCD), and their modeling requires assumptions about the nucleon wave 
function. 

When dealing with high-energy scattering where hadrons travel near the speed of light, 
the most natural tool to describe the nucleon is the light-cone Fock state expansion of 
its wave function. The hadronic state is decomposed in terms of iV-parton Fock states 
with coefficients representing the momentum light-cone wavefunction (LCWF) of the N 



partons 
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121 ] . In principle there is an infinite number of LCWFs in such an expansion. 



However, since the constituent quark models work so well phenomenologically, there must 
exist a light-cone description of the nucleon in which only the Fock components with a few 
partons are necessary. 

The quark distribution amplitudes defined in terms of hadron-to-vacuum transition ma- 
trix elements of non-local gauge-invariant light-cone operators and describing the three-quark 
component of the nucleon have been studied extensively in the literature [H 
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1 81 ] . In turn, the light-cone Fock expansion of a hadron state is completely defined by 



the matrix elements of a special class equal light-cone time quak-gluon operators between 



the QCD vacuum |0) and the hadron [19(. The authors of Refs. |20|, |2l|, |22J considered 



the wave-function amplitudes keeping full transverse-momentum dependence of partons and 
proposed a systematic way to enumerate independent amplitudes of a LCWF given a par- 
ticular parton combination. If one truncates the light-cone expansion of the proton state 
to the minimum Fock sector with just three valence quarks, one can write down the matrix 
elements of a class of three-quark light-cone operators which serve to define a complete set of 

rese matrix elements can be simplified using 
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22 1, and at the end one finds that six 



light-cone amplitudes within the truncation. T : 
color, flavor, spin and discrete symmetries {20 . 
amplitudes are needed to describe the three-quark sector of the nucleon LCWF. Depending 
on the imposed gauge fixing conditions such amplitudes are real or complex. In the latter 
case, the amplitudes contain final state interaction effects. 

With such amplitudes one can calculate nucleon observables. One could choose a phe- 
nomenological approach parametrizing and fitting them to data. Here, the wave-function 
amplitudes are modeled in a light-cone constituent quark model (CQM) which has been 
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successfully applied in the calculation of 
eralized parton distributions 
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;he electroweak properties of the nucleon 23j , gen 
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271 ] and spin densities 28| . This representation 



is well suited to disentangle the contribution from the different orbital angular momentum 
components of the nucleon wave function, and therefore to study the spin-spin and spin-orbit 
correlations encoded in the different TMDs. 

The paper is organized as follows. After a brief review of the definition of TMDs in 
Section [Tlland the resulting light-cone amplitudes when considering only three active quarks 
in Section [ITl] the proton wave function is constructed in Section HVl where the representation 
for the nucleon amplitudes in terms of the light-cone CQM is derived. Time-even TMDs 
are then calculated in Section |V] confining explicit expressions for the light-cone amplitude 
overlap representation of both time-even and time-odd TMDs in Appendix A. Numerical 
results are presented in Section IVIl and concluding remarks are given in Section IVHl 

II. TRANSVERSE MOMENTUM DEPENDENT PARTON DISTRIBUTIONS 



In this section we review the formalism for the definition of TMDs, following the conven- 
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tions of Refs. 
SIDIS is defined as 

$ij{x, k±, S) 



10]. The quark-quark distribution correlation function appearing in 



(2tt) 3 



where k + = xP + , and here and in the following we omit flavor indices. U is the Wilson link 



connecting the two quark fields and ensuring color gauge invariance of the correlator [29]. 
The target state is characterized by its four-momentum P and the covariant spin vector 
S = (S + , S~, S ± ) (P 2 = M 2 , S 2 = -1, P ■ S = 0), where 



p+ p- 

S + = A ^T> S- = -A—. 
M M 



(2) 



The TMDs enter the general decomposition of the correlator &ij(x, k±, S) which, at the 
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twist-two level, reads 

9lTlhV>+ + hi T - 75 

+ ih± [h ^\, (3) 

where n + and n_ are two lightlike vectors satisfying n + • n_ = 1, and = e and the 
transverse four- vectors are defined as vj_ = (0,0,v±). The nomenclature of the distribution 
functions follows closely that of Ref. [7|, sometimes referred to as "Amsterdam notation": 
/ refers to unpolarized target, g and h to longitudinally and transversely polarized target, 
respectively; a subscript 1 is given to the twist-two functions, subscripts L or T refer to the 
connection with the hadron spin being longitudinal or transverse and a symbol _L signals 



ght 

m 



the explicit presence of transverse momenta with an uncontracted index. Among the ei 

r 

distributions of Eq. (jHJ), the Boer- Mulders TMD hj; [8|] and the Sivers function 
are T-odd, i.e. they change sign under "naive time reversal", which is defined as usual 
time reversal, but without interchange of initial and final states. All the TMDs in Eq. ([3]) 
depend on x and k\. These functions can be individually isolated by performing traces of 
the correlator with suitable Dirac matrices. Using the abbreviation $t r l = Tr($r)/2, and 
restricting ourselves to the T-even TMDs, we have 



<^ + l(x,fc ± ) = h, 
*Vrt(x,k ± ) = Ag lL + ^^g rr . 



M 



M 



Ahi L + 



M 



h 



IT 



S ± h + A— h 1L + S ± — h 



IT- 



(4) 
(5) 

(6) 
(7) 



where j — x, y is a transverse index, and 



k 2 



(8) 
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The correlation function ^ 1+ ^(x,k±) is just the unpolarized quark distribution, which 
integrated over k^ gives the familiar light-cone momentum distribution fi(x). All the other 
TMDs characterize the strength of different spin-spin and spin-orbit correlations. The pre- 
cise form of this correlation is given by the prefactors of the TMDs in Eqs. (JS!)-©- In 
particular, the TMDs gn and hi describe the strength of a correlation between a longitu- 
dinal/transverse target polarization and a longitudinal (circular) /transverse (linear) parton 
polarization. After integration over kj_, they reduce to the helicity and transversity distri- 
butions, respectively. By definition the spin-orbit correlations described by giT, h^ L and h^ T 
involve the transverse parton momentum and the polarization of both the parton and the 
target, and vanish upon integration over k±. 

If one calculates these distributions in the light-cone gauge A + = using the advance 
boundary condition for the transverse com pon ent of the gauge field, the gauge links in the 

331. However, in this case the wave function 



quark-quark correlator can be ignored 
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amplitudes are not real, and apart from the structural information on the hadron it has also 
an imaginary phase mimicking the final state interactions. 

By using for the quark fields the canonical expansion in terms of Fock operators, one can 
write the T-even TMDs as 

mx,kl) = (PA\Y^q{(k)qx(k)\PA), (9) 

A 

Ag q 1L (x,kl) = (PA|^(-l) ( ^ A) gi(fc)gAW|PA), (10) 

A 

ik> ' S,) gl T (x,kD = (PS ± \J2(-^- X k{(k)q x Ck)\PSi_), (11) 

A 



M 



hl(x,kl) = \(PS x \qUk)qx(~k)\PS x ) +i(PSV\[^(\)]ql x (k)q x (k)\PSy)] , (12) 



2 

A 



A§ hf L \x, k\) = -(.y +1 (PA| ]r[sign(A)P+W%A(fc)|PA>, (13) 

A 

(kl - K) 



][2 h#(x, k 2 ± ) = [(PS x \qU~k)qx(k)\PS x } 

A 



-i(PSV\[sign(\)]ql x (k)q x (k)\PSy) , (14) 
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where x > 0, and q\(k) (q\(k)) is the annihilation (creation) operator of a quark of flavor q 
with helicity A and momentum k = (k + , k±). 

III. THREE-QUARK LIGHT-CONE AMPLITUDES 



In this section we reproduce the main results of Ref. |20j for the classification of the three- 
quark LCWF of the nucleon with helicity A in terms of the total parton light-cone helicity 
A, or, equivalently, in terms of the angular momentum projection l z — A — A which follows 
from angular momentum conservation. For a proton with helicity A = 1/2, the complete 
three-quark light-cone Fock expansion has the following structure 

The different angular-momentum components of the state in Eq. (1151) are given by 
|P T )' Z=0 = [d[l]d[2]d[3]U^(l,2,3) + zefk la k 2 ^ 2 \l,2,3) 



x^4(l) («} i (2)4 t (3) - 4 ; (2)4 T (3)) |0) , (16) 



|PT)^ =1 = y d[l]d[2]d[3] (^ (3) (1,2,3) + ^ (4) (1,2,3)) 

(^(1)^^(2)4,(3) - 4,(1)^,(2)4,(3)) |0> , (17) 

IP]) 1 *- 1 = J d[l]d[2]d[3](-)k~ ± ^ 5 \l,2,3) 

4(1) (^ T (2)4 T (3) - 4,(2)4,(3)) |0> , (18) 

\P |)^= 2 = / d[l]d[2]d[3] fci.*i.^ (B) (l,2,3) 
\/6 



X^4 4 (l) (4(2)4,(3) - ^(2)4,(3)) |0) , (19) 



where a, (3 = 1, 2 are transverse indexes and = kf ±/cf . In Eqs. ( IT6l) - (|T9l . the integration 
measures are defined as 

*ww - (i -|>) (|><) • < 20 ' 
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where Xi are the fraction of the longitudinal nucleon momentum carried by the quarks, 
and ki± are their transverse momenta. Furthermore, u\ x (ui\) and d\ x (di\) are creation 
(annihilation) operators of up and down quarks with helicity A and color i, respectively, 
and ipv) are functions of quark momenta with argument i representing Xi and fcjj_, and the 
dependence on the transverse momenta is of the form ki± ■ kj± only. 

The proton state with negative helicity is given in terms of the same wave function 
amplitudes except that the quark helicities are flipped, k x ± ik y become k x + ik y , and 
some signs are added 

\Pi) h=0 = J ' d[l]d[2]d[3] (-^«(l,2,3)+i^ , fc la fc^«(l,2,3)) 

x^^(l) (^(2)^(3) - 4(2)^(3)) |0) , (21) 

iPj)^- 1 = J d[l]d[2]d[3] (A; 1± ^(l,2,3) + A;-^ (4) (1,2,3)) 

X 7H K^Ht^t 3 ) ~ <(l)«t T (2)t4 T (3)) |0) , (22) 

\Pi) h=l = J d[l]d[2]d[3] (-)^ (5) (1,2,3) 

x^ilOO («5i(2)4i(3) - 4l( 2 X( 3 )) 1°) > (23) 

ipi)^=- 2 = /"d[i]d[2]d[3]H ^i^a^^) 



x^J T (l) (4 t (2)4 T (3) - ^} T (2)4 T (3)) |0> . (24) 

As the previous results are derived within light-front quantization, we implicitly assumed 
to work in the light-cone gauge A + = 0. However, this last condition does not fix the gauge 
completely, and additional boundary conditions must be specified. Depending on whether 
the additional gauge condition satifies time reversal or not, the wave function amplitudes 



are real or complex [20(] (see also [361 ]) 
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IV. NUCLEON WAVE FUNCTION IN A LIGHT-CONE CONSTITUENT 

QUARK MODEL 

In this section we derive the light-cone amplitudes ip^ in a light-cone CQM. Working 



in the so-called "uds" basis 



34 



351 ] the proton state is given in terms of a completely 



symmetrized wave function of the form 



|Pt) = \P Vuud + \P Vudu + \P t) 



duu 



(25) 



In this symmetrization, the state \P ]) u du is obtained from \P ]) UU( i by interchanging the 
second and third spin and space coordinates as well as the indicated quark type, with a 
similar interchange of the first and third coordinates for \P ])duu- 

Following the derivation outlined in Ref. 24J, we find that the uud component of the 
light-cone state of the proton can be written as 



/ ijk 
ci[lM2M3]*^^ (26) 



where assuming SU(6) spin-flavor symmetry the LCWF ^tud(i x ^ ^JL<; \}) ls given by 



^tu f d\{^ku;X l }) = ^({x i ,k 1 _ i })^ A (X 1 ,X 2 ,X 3 ). 



1 = 



(27) 



In Eq. f l27|) the momentum dependent wave function is defined as 

nl/2 



^({x h k ±i }) = 2(2vr) 



1 U1U2UJ3 



i/j({xi,k ±i }), 



(28) 



_M X 1 X 2 X 3 

with ip({xi,kj_i}) symmetric under exchange of the momenta of any quark pairs, u>i the 
free-quark energy, and M = Y^=i u i the mass of the non-interacting three-quark system. 
The spin dependent part in Eq. ( |27l) is given by 



$a(Ai,A 2 ,A 3 ) = (l/2,/Ui;l/2,/i 2 |l,/i 1 + / u 2 )(l, / u 1 + / u 2 ;l/2,/i 3 |l/2,A) 

xDlitiRcM, k L1 ))D]! 2 2 * 2 {R cf {x 2 , k^D'JXiRcfixs, fcjLs)), (29) 



where D 1 x ^(R c f(x, k±)) is a matrix element of the Melosh rotation R c f 37], 

D^ 2 (R cf (x,k ± )) = (X\R cf (x,k ± )\fx) 

. , m + xM — icr ■ (z X k±) 

\X\ / — \w- 

m + xM ) 2 + k]_ 



(30) 
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The Melosh rotation corresponds to the unitary transformation which converts the Pauli 
spinors of the quark in the nucleon rest-frame to the light-front spinor. In particular, the 
spin wave function of Eq. (1291) is obtained from the transformation of the non-relativistic 
spin wave function with zero orbital angular momentum component. The relativistic spin 
effects are immediately evident in the presence of the spin-flip term icr ■ (z x k±) in Eq. fl30l ). 
Such a term generates non-zero orbital angular momentum, and, as a consequence of total 
angular momentum conservation, total quark helicity different from the nucleon helicity. 
Making explicit the dependence on the quark helicities, the spin wave function of Eq. (129]) 
takes the following values: 

$t(T,U) 
MU.T) 
<MU,t) 
<MT,U) 
<MU,I) 
<MU,T) 
MU,T) 
<MU,I) 



n 

i 


1 




1 








n 
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1 




1 


y/N( Xi , 






n 
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1 




1 


y/N{ Xi , 




V6 


n 

i 


1 




1 




k±i) 


V6 


n 

i 


1 




1 




k±i) 


V6 


n 

i 


1 




1 




k±i) 


V6 


n 

i 


1 




1 




k±i) 


V6 


n 

i 


1 




1 


y/N(x i} 


k±i) 


V6 



(2a!a 2 a3 + aifcjfcj" + c^fcjffcj - ), 


(31) 


(— a!a 2 a 3 + a^k^k 2 — 2ai/cjA;j), 


(32) 


(— axa 2 a3 + a^kfk^ — 2a 2 kfk^), 


(33) 


(aia 2 k^ — k^k 2 k^ — 2a\a^k^\ 


(34) 


(—k^k^k^ + 0102^ — 2a 2 a^,k^), 


(35) 


(a 2 a^k^ + aio^/c^ + 2k^k 2 k^), 


(36) 


(— a 1 a 3 A; 2 " — 0203/^ + 20102^3"), 


(37) 


(—a 2 kfk^ — a\k 2 k^ + 2a^kfk 2 ), 


(38) 



where a« = (m + XjM ), and N(xi, k±{) = {{m + XiM ) 2 + k\ { ]. 

Taking into account the quark-helicity dependence in Eqs. (l3~i~!) - (!38"l) . the nucleon state 
can be mapped out into the different angular momentum components of Eq. (fTSl) . After 
straightforward algebra, one finds the following representation for the nucleon amplitudes 
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in the light-cone CQM 

^ (1) (1,2,3) = 4>({x h k ±l }) 

x TT = — (-aia 2 a 3 + a 3 fci ± • fc 2± + 2aifcn. • fc 2 ± + 2ai&2 ± ), 

7V%,y V3 

(39) 

V (2) (l,2,3) = ^{x^kuVH 1 J^(a 3 + 2a 1 ), (40) 

, y/N(Xi, k_Li) V3 

^(1,2,3) = -^(fo.fc^JI 1 -Ua^ + fc^), (41) 

V> (4) (1,2,3) = -^({gi,fc±i»TT , / , ^ ( a i a 2 + 2Q 3 Q i- fc i±- 2fc i±- fc 2±), (42) 

j VN{x t ,k ±i ) V3 

^(1,2,3) = ^(K fcii}) J] TwT i , 4fK Q 3), (43) 

1 \/N(xi,k±i) V3 

^ (6) (1,2,3) = ^({x t ,fc ±t })n 1 = ^a 2 . (44) 

2 y/N(Xi, k_Li) V3 

The results in Eqs. (139 1) - 0441) follow from the spin and orbital angular momentum structure 
generated from the Melosh rotations, and are independent on the functional form of the 
momentum dependent wave function which we assumed symmetric under permutation of 
any quark pairs. 



V. TMDS IN A LIGHT-CONE CONSTITUENT QUARK MODEL 

General expressions of the TMDs can be derived under the assumption that the light- 
cone Fock expansion can be truncated to just the three valence quark contribution so that 
the nucleon wave function can be expressed in terms of six amplitudes. In terms of matrix 
element between proton states with different orbital angular momentum components, TMDs 
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are given by the following expressions 

A A 

A A 

g q iL {X,kl) = h=0 (P t \Og\P |) i2=0 + h = l (P]\O g \P]) h = l 

+ l * = - X (P t \Og\P ]) lZ= ~ l + h=2 (P T \Og\P }) lZ=2 ) (46) 

hl(x, k\) = Re[<*=°(P | |g|g T |P T>^ = °] + 2Re[^=- 1 (P | l^tl^ i>'* =_1 ], ( 4 ?) 

4%, k\) = -Re[^ =1 (P T |g|g T l^ i) 1 ^) ~ 2Re[<*=°(P T |g{? T l^ i)^ 2 ], (48) 
*i) = ™ (Re[ /z=1 (P T |g{g T |P T)' z=0 ] - Flm^P T |g|g T l^ T)' z= °] 

+rRe[' z=0 (P t |g{g T l p T) iz=_1 ] - k y Im[ l * =0 (P T |gJg T |P T)' z=_1 ] 

+PRe[ , '= 2 (P T \q\q<t\P T) iz=1 ] - ^Im[^= 2 (P | |<?{g T |P T) iz=1 ]) , (49) 
«it(*, fc D = (^Re[^=°(P T |0,|P l)^ 1 ] + FIm[^=°(P T \O a \P i) lz= ~ 1 } 

+FRe[' z= - 2 (P | \O g \P ]) lz= ' 1 } + FIm[^=- 2 (P | \Og\P ]) h= ' 1 ]){^) 



where = Ea(-1) 1/2_A <7a<?a. 

The unpolarized TMD /f, the helicity TMD g^ L , and the transversity TMD h\ in 
Eqs. (1451)- (1471) involve matrix elements which are all diagonal in the orbital angular mo- 
mentum, but probe different transverse momentum and helicity correlations of the quarks 
inside in the nucleon. In particular, ff is defined in terms of the momentum density opera- 
tor, g\ L is sensitive to the difference of right and left quark-helicity. Viceversa, h\ involves a 
chiral-odd operator with a quark-helicity flip compensated by a flip of the nucleon helicity 
in the same direction. The same chiral-odd operator enters the definition of h^j, and h\l in 
Eqs. (148 p and (j4T)l) . respectively. In the case of the nucleon helicity flips in the direction 
opposite to the quark helicity, with a mismatch of the orbital angular momentum between 
the initial and final nucleon state of Al z = 2, whereas hfl is diagonal in the nucleon helicity, 
with the quark-helicity flip inducing a change by one unit in the orbital angular momentum 
of the initial and final nucleon state. Finally, g\ T is defined in terms of the same helicity 
operator which enters the definition of gf L , but this time the nucleon helicity flips, with a 
transfer of orbital angular momentum by one unit. 

If one inserts in Eqs. (I45l) - fl50l) the three-quark light-cone amplitudes introduced in 
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Eqs. (17)-(19) and (22)-(25) of Section HTT1 one obtains the general expressions of the TMDs 
collected in Appendix A. These formulas can be worked out in terms of the explicit repre- 
sentation of the light-cone amplitudes obtained in the light-cone CQM of Section HVl As 
we neglected gluon degrees of freedom, the amplitudes in the light-cone CQM are pure real 
functions and lead to the following results for the T-even TMDs 



ft{x,kl) = N«f d[l]d[2]d[3] > /xIii^<5(A;-A;3)|V'({a:i},{fc± > i})| 2 , 

/( rn i r ]\/f \2 _ t.2 
d[i}d[2}d[3}^E^- 3 s(k - k 3 m{ Xt }, {fc^Di 2 i J:, A , 
(to + xMq) 2 + kj 

gl T {x,kl) = P* f d[l]d[2]d[3]V5[i^(fc ~ {k±M % 2M }™^ X 2 M °L , 

J [rn + xiVIq) -\- Kj_ 



(51) 



(52) 



h\{x, k 



P" 



d[l]d[2}d[3}VwT 3 6(k - kM({ Xi }, {fc ± ,})| 2 (m \ X , M " ? l2 

[rn + xM ) 2 + fcj 



(53) 



.pi 



d[l]d[2]d[3]^/xTx^r 3 5(k - k 3 )\^({x t }, {k Xti })\< 



(54) 



2M 2 



(m + xMof + kY 



(55) 



hf«(x,kl) = -Pi I d[l]d[2]d[3]Vx^T 3 5(k - k 3 )m{ Xl }, {k ± ,})\ % 2 ^ (m A# + ^2 . 

J (TO + xMq) z + fe_L 

(56) 

4 1 

where 5(k — k 3 ) = 5(x — x 3 )5(k± — k± 3 ), and the flavor factors P u = -, P d = — , N u = 2 

3 3 

and N d = 1 are dictated by SU(6) symmetry. 

By inspection the above TMDs satisfy the following relations: 



2hl(x,kl) = g q 1L {x,kl) + ^- q n{x,k\ 



P q o 2 &i 2 

■jy-fi(x,k ± ) = h^XykjJ — ^j^h 1T (x,k ± ), 
h^ix^k^ = —g^ T [x,k ± ). 



(57) 

(58) 
(59) 



Eq. (j57p is a generalization of analogous relations discussed in Ref. |26l . |38| | and was also 
rederived together with Eq. (|58|) in Ref. [39J. Eq. (|59|) was already found in the diquark 
spectator model of Ref. 48J. In QCD TMDs should be all independent of each other. The 
limitation to three valence quarks implies that out of the six TMDs f\, g±L, gir, hi, Jv^r , 
h^ L only three are linearly independent. A similar situation occurs with the bag model 39]. 



13 



In the bag model there are only S- and P-wave components of the proton wave function, 
whereas here also a D-wave contributes. However, the relations fl57|) - fl59|) do not depend on 
the different components of orbital angular momentum. Their specific form is a consequence 
of the imposed SU(6) symmetry which allows us to factorize the momentum dependent wave 
function from the effects of the Melosh rotation acting in the spin space and producing 
different factors for the different TMDs, Eqs. (loTj) - f[o'6"j) . In the diquark spectator model 
the relations (15Tj) and f)58p hold only for the separate scalar and axial contributions, while 
Eq. (159|) is verified more generally for both u and d flavors. 

Concerning what has been called the 'pretzelosity' distribution, h^ T |39j, from Eqs. fl5Tl) . 
( 152|) . and (1551) one easily verifies the positivity condition 49 1 



2M 2 h lT 



(60) 



Furthermore, subtracting the relations (|57|) . (1581) from each other yields 

k 2 



h ±q (T k z 



gf L (x,ki)-hl(x,ki). 



(61) 



This result was already found in Ref . 39J . Integrating out transverse momenta and going to 
the non-relativistic limit where helicity and transversity distributions coincide, one finds that 
the first moment of h^ T vanish identically. Thus, relation ( ISTT) supports the statement that 
hj; T is a measure of relativistic effects. Relativity, responsible for a chiral-odd transversity 
distribution differing from a chiral-even helicity distribution, exhibits the chirally odd nature 
of hi T . This is confirmed by the following relation that is also satisfied within our model: 

3 



h 



IT 1 



X) 



(1-x) 



;H q T (x,0,0), 



(62) 



where the transverse moments of are defined as 



X 



I MM 



and Hj,(x, 0, 0) is the forward limit of a chiral-odd generalized parton distri 



(63) 



bution (GPD) 



2fl 



4l|). Eq. (jn2D 



occurring in the case of parton and nucleon helicity flip (see, e.g., Refs. 
was first found in Ref. 40J] to hold for the scalar diquark model and in a quark target model 
of the nucleon as a particular case of a general relation between the moments of TMDs and 
the moments of GPDs. 
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By integrating over x the first moment of hj- T , from Eq. (I6T1) one obtains 

/ dx h^ ±q (x) = Aq- Sq, (64) 

where Aq and Sq are the axial and tensor charges that measure, for each flavor q, the net 
number of longitudinally polarized valence quarks in a longitudinally polarized nucleon and 
the net number of transversely polarized valence quarks in a transversely polarized nucleon, 
respectively. 



VI. RESULTS 



48j (see also [50j, |5l| ) using 



The full list of T-even quark TMDs was computed in Ref. 
the diquark spectator model with scalar and axial-vector diquark. The analytic form of 
TMDs was also derived in a quark target model 40J]. Here, the formalism described in the 
previous sections is applied in the following to a specific CQM adopting a power-law form 
for the momentum dependent part of the light-cone wave function 42], i.e. 

1/2 



^({x i} k ±i }) = 2(2tt) 



1 UJ1U2UJ3 



N' 



(M 2 + /3 2 )^ 



(65) 



_M x 1 x 2 x 3 

with N' a normalization factor. In Eq. (165|) . the scale (3, the parameter 7 for the power-law 
behaviour, and the quark mass m are taken from Ref . [421] . i.e. (3 = 0.607 GeV, 7 = 3.4 
and m = 0.267 GeV. According to the analysis of Ref. 43J these values lead to a very good 
description of many baryonic properties. 

The results for ff, g\ L and h\ are shown in Fig.[TJ They are consistent with those obtained 
in Ref. Q for the corresponding PDFs indicating that ft and ff have the same (positive) 
sign, whereas g\ L and h\ have opposite sign with respect to gf L and hf, respectively. In 
addition, the size of the TMDs for d quarks is smaller than that for u quarks according to 
the flavor dependence of TMDs through the factor P q in Eqs. ( f52l) -([56l). It is remarkable 
that the k 2 ± dependence cannot be factorized as a Gaussian function as often assumed. In 
all distributions the maximum (minimum) is around x = 0.2 at k\ = and moves to higher 
values of x with increasing k\. 

The TMDs hf^ and hfj! are shown in Fig. [2j The size of both and h^f is much larger 
than that of ff, g\ L and h\, a result in qualitative agreement with that obtained in the bag 
model 40]. In particular, comparing h\ and one deduces that the quark helicity flip 
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is more favored in the case of a longitudinally polarized nucleon with a transfer of orbital 
angular momentum between initial and final states. The shapes of the ^-distributions of 
h\ and h x £ are similar, but with opposite sign. On the other hand, hyp has a narrower 
x-distribution with a faster fall-off in k 2 ± . The pretzelosity hfp contributes when the quark 
and nucleon helicity flip in opposite directions. It then requires an overlap between wave 
function components that differ by two units of orbital angular momentum, either a PP or an 
SD interference (see Eq. ( I48p ). The different partial wave contributions to hfp are plotted in 
Fig. [3] where one may notice the importance of also including the D-wave component which 
is absent, e.g., in the bag model and the diquark spectator model. While in the case of u 
quarks the PP and SD interference terms add with the same sign, in the case of d quarks 
they have opposite sign, indicating that the SU(6) relation between u and d contributions, 
hj?p = —4hyp, is valid for the total result but not for the partial wave contributions. 

The results presented in Figs. [J and [2] are qualitatively similar also to those obtained with 



the diquark spectator model 



48] . With respect to the results obtained in Ref . [48] with scalar 



and axial masses M R = 0.6 and 0.8 GeV, respectively, and with a cut-off A = 0.5 GeV, TMDs 
calculated with the light-cone CQM are peaked at smaller values of x with broader x and 
k\ distributions. As a consequence, the transverse moments of the pretzelosity distribution 
may be rather different in different models, as can be seen in Fig. H] where the light-cone 
CQM results with the momentum wave function (1651) are compared with similar calculations 



with another momentum wave function derived within the hypercentral model 46|, |47j and 



with results of the bag model and the spectator model of Ref. [48]. This sensitivity to the 
adopted model suggests that new data could give useful insights to model the momentum 
dependence of the nucleon wave function. 

According to Eq. (}59l) . g\ T is just the opposite of h x £. Both TMDS g\ T and h x £ involve 
matrix elements between states that differ by one unit of angular momentum and one finds 
here h\p < < h\£ and g XT > > gf T , in agreement with some expectation 44] . On 
the contrary, in our model h yp and hy/ are obtained with the reversed sign predicted by 
qualitative arguments in Ref. 44] but in agreement with the result within the bag model 40 ] 
and the diquark spectator model 



48]. 

In Ref. [45] approximate relations among TMDs were studied. Taking advantage of 
the QCD equation of motion and neglecting twist-three TMDs, Wandzura-Wilczek-type 
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approximations were proposed for the transverse moments of g\ T and i.e. 

ml dy 



g$\x) * x ! iLgXy), (66) 

j x y 

h^\x) « - x * f d JL h \{y). (67) 



According to Eq. (1551) . and fojx should be equal and with opposite sign in the 

present model, whereas the helicity and transversity distributions, and h\ respectively, 
are rather different 



381 ] . As a consequence, in Fig. [5jone may appreciate how good is the 
Wandzura-Wilczek-type approximation when considering an SU(6) symmetric model as the 
one adopted here. In fact, the approximation of neglecting twist-three contributions works 
better for hf^ q than for g^ q '. In any case, the model results support that the estimates 



for spin observables in SIDIS made in Refs. 45|, |52|] on the basis of the approximations (66) 
and (67) have a useful accuracy. 

The presence of a significant orbital angular momentum component is suggesting a non- 



spherical shape of the nucleon. According to Ref. [53( there is an infinite variety of obtainable 
shapes depending on the contribution of the pretzelosity h^ T . They can be found by look- 
ing at a suitable spin- dependent quark density, Prelt, in a nucleon state polarized in the 
transverse direction St either parallel or antiparallel to a given direction n. The transverse 
shapes of the nucleon are then derived from the following relation: 

p REL r(fc±,n)/M h x (kl) k 2 ± Mr( k l) fRa ^ 

T71Y\ = 1 + j / i 2 \ cos <p n -\- — — cos(20 - n ) , 68 

fx{k\) h{k\) 2M 2 f^kl) 

where <fi is the angle between k± and St and n is the angle between n and St- A tilde is 
placed over a given quantity to define the ^-integrated result, e.g., 

h(kl) = j dxh{x,kl). (69) 

Assuming a struck u quark, the transverse shapes of the proton are shown in Fig. [6] for 
St parallel to n, <p n = 0, and in Fig. [7] for St antiparallel to n, <p n = ir. The corresponding 
results assuming a struck d quark are shown in Figs. [8] and [9j respectively. In our model /" 
(fi) and (hf) are of the same (opposite) sign and similar size, so that the contribution 
of the first two terms on the rhs of Eq. (|68l) tend to cancel each other for <ft n = it (<f) n = 0) 
emphasizing the role of the pretzelosity in producing deformation. For u (d) quarks the last 
term in Eq. (IBS"]) is negative (positive) for (ft = 4> n = and its size increases (reduces) with 
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the inclusion of the D-wave. This explains the larger transverse deformation in the direction 
antiparallel (parallel) to St for a struck u (d) quark, with a more significant effect in the 
case of the u quark. 



VII. CONCLUDING REMARKS 



Quite general expressions have been derived for the transverse momentum dependent 
parton distributions in a light-cone description of the nucleon where the Fock expansion 
is truncated to consider only valence quarks. They are given in terms of matrix elements 
between nucleon states with different orbital angular momentum so that one can immediately 
appreciate the origin of each individual TMD. In particular, combining the total parton light- 
cone helicity with the nucleon helicity, angular momentum conservation requires the angular 
momentum projection of the nucleon to run from —1 to 2, thus imposing S-, P- and D-wave 
components in the nucleon wave function. 

The complete three-quark light-cone wave function involves six amplitudes that can be 
either real or complex depending on the gauge fixing conditions 20Ll2l|. When complex, they 
incorporate the effects of final state interactions. The light-cone amplitudes have been used 
to obtain a model independent light-cone amplitude overlap representation of the T-even and 
T-odd TMDs, which emphasizes the role of the different angular momentum components. 
In the present approach where gluons are not included, the six amplitudes are real and 
have been constructed by assuming a light-cone constituent quark model with SU(6) spin- 
flavor symmetry. Analytic formulae have been derived for the T-even TMDs in terms of 
the momentum dependent part of the LCWF. By simple inspection of such formulae we 
have found that among the six twist-two T-even TMDs there are three relations, so that 
only three TMDs are in fact independent when assuming SU(6) spin-flavor symmetry in 
the three-quark sector of the Fock expansion. Two of such relations, Eqs. (1571) and fl58l) . 
were already known 38|, |39(, whereas the third one, Eq. (1591) . was first found in the diquark 
spectator model 48( and shown here to be valid in a larger class of relativistic quark models 
with SU(6) symmetry. 

Numerical calculations of the T-even TMDs have been presented by adopting a power- 
law form of the momentum dependent part of the LCWF. All distributions have either a 
maximum or minimum at x ~ 0.2 and k\ = moving to higher values of x with increasing 
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k\. This indicates that the usual Gaussian ansatz for a factorized k\ dependence of TMDs is 
not adequate. In any case, the formalism described here is well suited to study the effects of 
the momentum dependence of the nucleon wave function, in particular in phenomenological 
applications to SIDIS and DY processes. 

The role of the different orbital angular momentum components in the nucleon wave 
function is best appreciated by looking at the so-called pretzelosity, h^ T , and its effect on 
the spin- dependent quark density producing a nonspherical shape of the nucleon [53]- A 
larger deformation in the direction antiparallel (parallel) to the transverse spin St is found 
for a struck u (d) quark, with a significant sensitivity to the presence of a D-wave component 
in the nucleon wave function. 
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Appendix A: LCWF OVERLAP REPRESENTATION OF TMDS 

Explicit expressions of the T-even TMDs in terms of overlap of the light-cone amplitudes 
of Section II III are given in this Appendix. The results for g\ T and h x £ have already been 
derived in Ref. [2(3], while the results for the other T-even TMDs are given here for the first 
time. However, our expressions for g\ T and h\l differ from the previous results of Ref. 
because we have corrected the wave function component with tp^ according to Ref. 21 1 

The TMD ff(x, k\) reads 



•,f. 2(J 



f!{x,k 2 ± ) = J d{l]d[2)d[3}^^ F", (Al) 
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where for the up quark we have 

T u = 2 5 3 (k - fci)^ (1>2) (l, 2, 3)$ (1 ' 2) *(1, 2, 3) + ^ (1 ' 2) *(3, 2, 1)] 
+ [5 3 (A; - fc 3 ) + 5 3 {k - k 2 )} [^ 2 \l, 2, 3)^*(1, 2, 3)] 
+2 <5 3 (A; - fc 2 )V> (3 ' 4) (l, 2, 3)[^ {3,4) *(1, 2, 3) + ^ 3 > 4 >*(1, 3, 2)] 
+ [5 3 (k - k,) + 5 3 (k- k 3 )} [^(l, 2, 3ty< 3 ' 4 >(l, 2, 3)] 

+[S(k - h) + 6(k - k 2 )][ft 5 \l, 2, 3) + ^(2, 1, 3)] 0< B >(1, 2, 3) - # 8 >'(1, 3, 2)] 
+[8( k - h) + 6(k - k 3 )][fi*\l, 2, 3) + ^ {5) (3, 2, 1)] 0< B >*(1, 2, 3) - ^< 5 >'(1, 3, 2)] 

+ [6 ^ k _ fcl ) + j3 (fc _ fcs)] ^(6) (1> 2> 3) + ^(6) (3? 2> 1)][ ^(6)* (1; 2) 3) _ ^(6> (1> 3 2)]; 

+ [5 3 (A; - fcO + 5 3 (k - k 2 )\ 2, 3) + ^ (6) (2, 1, 3)]0< fl >'(l, 2, 3) - # 8 >*(1, 3, 2)],(A2) 

with S 3 (k — ki) = S(x — Xi)S(k± — k i± ) and 

^ 2 )(1, 2, 3) = 2, 3) + iefk la k 2p ^ 2 \l, 2, 3), (A3) 

<P' 4 )(1, 2, 3) = fc+^l, 2, 3) + kl^\l, 2, 3), (A4) 

^ 5 )(1,2,3) = A; 2 -^ 5 )(1,2,3), (A5) 

^ 8 )(l,2,3) = fcffc^ 6 )(l,2,3). (A6) 

For the down quark, we obtain 

T d = 5 3 (k - fc 2 )^ (1>2) (l, 2, 3)$ (1 > 2) *(1, 2, 3) + ^ (1 ' 2) *(3, 2, 1)] 
+ < 5 s (fc-fc 3 )[^ (1 ' 2) (l,2,3)^ 1 ' 2 )*(l,2,3)] 
+5 3 (k - 2, 3)[^ 3 - 4 )*(l, 2, 3) + <p- 4 )*(l, 3, 2)] 

+5 3 (A;-A;3)[^( 3 ' 4 )(1,2,3)^ 3 ' 4 )*(1,2,3)] 

+5 3 (A; - A; 3 )[^ (5) (l, 2, 3) + ^ {5) (2, 1, 3)] $< B >*(1, 2, 3) - ^ (5) *(1, 3, 2)] 
+5 3 (A; - fc 2 )0< 5 >(l, 2, 3) + ^ {5) (3, 2, 1)] 0< 5 >*(1, 2, 3) - ^>(1, 3, 2)] 
+5 3 (A; - k 2 )[ft 6 \l, 2, 3) + ^(3, 2, 1)] 0< 6 >*(1, 2, 3) - ^ 6 )*(1, 3, 2)], 
+5 3 (k - k 3 )[^(l, 2, 3) + ^ (6) (2, 1, 3)] 0<«>*(1, 2, 3) - # 8 >'(1, 3, 2)]. 

The result for gf L (x, k 2 ± ) reads 

g q 1L (x,k 2 ± ) = f d[l]d[2]d[3]y/ww Gl 



(A7) 
(A8) 
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where the function Q\ for the up quark is 



Ql = 2 5 3 (k - h)^ 2 \l, 2, 3)[^ (1 ' 2) *(1, 2, 3) + ^*(3, 2, 1)] 
+ [5 3 (k - fc s ) - 5 3 (k - k 2 )][^ 2 \l, 2, 3)^ 2 >(1, 2, 3)] 
-2 5 3 (A; - A; 2 )^ {3 ' 4) (1, 2, 3)$< 3 ' 4 >(1, 2, 3) + ^ 3 - 4 >*(l, 3, 2)] 
+ {5 3 (k - h) - 5 3 (k - A; 3 )][^ (3 ' 4) (l, 2, 3)^ 3 ' 4 >(1, 2, 3)] 

+ [S{k _ fcl ) + _ fc 2 )] [^(5) (1) 2 , 3) + ^ (5) (2, 1, 3)] $< 5 >*(1, 2, 3) - ^ 5 >(1, 3, 2)] 
+ [S{k _ + ^ _ fca )] [^(5) (1; 2 , 3) + ^ (5) (3, 2, 1)] 0< 5 >*(1, 2, 3) - # B >*(1, 3, 2)] 
-[5 3 (k - h) + 5 3 (k - A; 3 )][^ (6) (l, 2, 3) + ^(3, 2, 1)] [ft 6 >(l, 2, 3) - ^< 6 >*(1, 3, 2)], 
_[ 5 3 (A; _ + 5 3 {k _ fc2 )][^(«)(i > 2 , 3) + ^(2, 1, 3)] 0< 6 >*(1, 2, 3) - ^< 6 >'(1, 3, 2)],(A9) 



G d L = S 3 (k - k 2 )^ 2 \l, 2, 3)[^ 2 >(1, 2, 3) + ^< 1 > 2 >(3, 2, 1)] 

+ < 5 3 (fc-fc 3 )[^ (1 ' 2) (l,2,3)^ 1 ' 2 )*(l > 2,3)] 

+5 3 (k - h)^ 4 \l, 2, 3)[^ (3 ' 4) *(1, 2, 3) + ^ 2 )*(1, 3, 2)] 

-5 3 (k-k 3 )[^ 4 \l,2,3)^ 4 >(l,2,3)} 

+5 3 (k - k 3 ) 2, 3) + ^ (5) (2, 1, 3)][^ (5) *(1, 2, 3) - ^ {5) *(1, 3, 2)] 

+5 3 (A; - k 2 ) $< 6 >(1, 2, 3) + ^(3, 2, 1)][^ (5) *(1, 2, 3) - ^ 6 >*(1, 3, 2)] 

-5 3 (k - k 2 )[^\l, 2, 3) + ^(3, 2, 1)]0< 6 >'(1, 2, 3) - ^>(1, 3, 2)], 

-5 3 (k - h) (1, 2, 3) + (2, 1, 3)] [# 8 >*(1, 2, 3) - ^*(1, 3, 2)]. (A10) 



and for the down quark is 



For h\ we have 




(AH) 



where for the up quark 



n 



2 {s 3 (k - h)[^ 2 \l, 3, 2) + ^> 2 )(2, 3, 1)]^*(1, 2, 3) 
-5(k - h)[i^(l, 2, 3) + ^(2, 1, 3)]^ 3 > 4 >(1, 2, 3) 
+5 3 (k - fca) [^(1, 3, 2) + ^(2, 3, 1)]^ 3 ' 4 )*(2, 1,3)}, 



(A12) 



with 



^(l, 2, 3) = ^(1,2,3) - 
^(1,2,3) = A; 2 V 5) (1 ) 2,3). 



iefk la k 2 ^ 2 \l,2,3), 



(A13) 
(A14) 
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For the down quark, we obtain 

H d = {-5 3 (k- fc 3 )^ 1 ' 2 )(l,2,3)^^*(2,l,3) 

+25(k - fc s ) 2, 3) + ^(2, 1, 3)]^ 3 < 4 )*(3, 2, 1) 

-25 3 (A; - k 2 ) [^(1, 2, 3) + <F(3, 2, l)$< 3 - 4 >(2, 1, 3) } 
= 5\k - k 3 ) Re {-^ (1 ' 2) (1, 2, 3)^*(2, 1, 3) 
+2[^(1, 2, 3) + ^(2, 1, 3)]^ (3,4) *(3, 2, 1) 

-2[j®(l, 3, 2) + ^(2, 3, 1)]^( 3 ' 4 >(3, 1,2)}. (A15) 
The result for g\ T is 

r 

g q 1T (x, k\) = £r y rf[l]rf[2]rf[3]v^i^Re^. (A16) 

The function ^ for up quarks is 

0? = 2{<5 3 (A;-A; 1 )# 2) *(1,2,3)^ 3 ' 4) (2,1,3) 

-<5 3 (A; - fc 2 )^ 1,2) (l, 2, 3)^ {3 ' 4) *(2, 1, 3) 

+ [5 z {k _ fcl ) + S *(k _ fc 3 )]^ 2) *(l, 2, 3)[</f 4) (2, 1, 3) + </f 4) (2, 3, 1)] 
+ [5%k - h) + 5 3 (k - A; 2 )][^ (5) *(l, 2, 3)^ (6+) (2, 1, 3) 

+^ (5) *(1,3,2)^ 6+ ')(3,2,1)]}, (A17) 

where 

Vf' 2) (l,2,3) = ^(1,2,3)^-^(^-^(1,2,3), (A18) 
^f' 4) (l,2,3) = fcj^(l,2,3)+^ 4 )(l,2,3), (A19) 
^ 6+ ')(l,2,3) = (k 2 2± k 3± -k ± + k 2 3± k 2± -k ± ) (V> (6) (2,l,3)+V (6) (3,l,2)) 

+fc 2 ± ^2± • fc ± ^ (6) (2, 3, 1) + k 2 1± k 3± • fc ± V (6) (3, 2, 1), (A20) 
^ {6+) (l,2,3) = -k 2 1± k 3± -k ± ^ 6 \l,2,3) 

+ (k 2 2± k 3± ■ k ± + fc 2 ± fc 2± • k ± ) ^ (6) (2, 1, 3) 

+k 2 1± k 2± ■ k± (^ (6) (1, 3, 2) + ^ 6 >(2, 3, 1)) , (A21) 
with i,j = x,y. For the down quark we obtain 

G d T = 2 {^(A; - fc 3 )^ (1 ' 2) *(l, 2, 3)^ (3 ' 4) (2, 1, 3) 

- 5\k - fc 2 )^ 1,2) (l, 2, 3) [^ (3 ' 4) *(2, 1, 3) + ^ (3 ' 4) *(2, 3, 1)] 

+5 3 (k - k 3 ) [^ 5 >(1, 2, 3)^ 6+ )(2, 1, 3) + V (5) *(l, 3, 2)^ 6 +')(3, 2, 1)]} . (A22) 
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As outlined in Ref. 20j the imaginary part of the function Qj, in Eq. 



sponding LCWF overlap representation for the T-odd Sivers TMD f^? 
The result for hf q is 



HI) gives the corre- 



30 . 



h 1 q (x, k]_) 



M 



Ak 



, , d[l]d[2]d[3]y/x 1 x 2 x 3 ReHi q , 



where the function H^ q for the up quark is 



5\k - fcOV^'Cl, 2, 3) [Vf*'*(3, 2, 1) + Vf* '*(3, 1, 2)] 
-5 3 (fc - fc a ) [^ 1 ' 2) '(1, 3, 2) + ^ 1>2) '(2, 3, l)]Vf 4) *(1, 2, 3) 
-5 3 (A; - fc 2 ) [V^CL, 2, 3) + ^(2, 1, 3)]^ 1,8) *(1, 2, 3) 
+5 3 (fc - ki) [^ {6 ' +) *(l, 2, 3)^(1, 2, 3) + ^ 6 -+>*(2, 1, 3)V (4) (1, 2, 3)], 



(3,4)*/ 



(3,4)*, 



(A23) 



(A24) 



where 



^(W) = fc^ (5) (l,2,3)- ^(1,3,2), 
# 2) '(1,2,3) = F'^ (1) (l,2,3) + e^F(A;^|-^^)^( 2 )(l,2,3). 



(A25) 
(A26) 



The result for the down quark is 



Ui d = 2 \6 3 (k - Wj 1,2) '(l, 2, 3)V>f 4) *(1, 2, 3) 



-5 3 (A; - fc 2 )Vf 2) *(1, 2, 3) [^(1, 2, 3) + ^(3, 2, 1)] 

+S 3 (k - fci)[^ 6 ' + >'*(l, 2, 3)V (3) (1, 2, 3) - ^ {6 ' +) *(2, 3, 1)^ (4) (1, 2, 3)]} . (A27) 

Note that by taking the imaginary part of the function 7i^ 9 in Eq. ( 1A23I) one obtains the 
corresponding results for the T-odd Boer- Mulders function hi q 
Finally, for the TMD we obtain 



hyp(x, k ± ) 



2M 2 

U2 _ U2 



d[l]d\2]d\5\y/x\X2Xz ReH^ q , 



(A28) 



where TL^ q is given by 



n 



±u 



2 \5(k - fc 2 )^ (3 ' 4) (l, 2, 3)[^ (3 ' 4) *(3, 2, 1) + ^ 3 ' 4 >(3, 1, 2)] 



-5 3 (k - h) 2, 3) + ^ 6 >(3, 2, 1)]^>*(3, 1, 2) 

-<5 3 (A; - k 2 )[^\l, 2, 3) + ^ 6 >(2, 1, 3)]^*(1, 2, 3) 



(A29) 
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with 

^(1, 2, 3) = k^\l, 2, 3) + fc^l, 2, 3). (A30) 
For the down quark, one finds: 



H^ d = - {8 3 (k - k 3 )[ft 3 ' 4 \l,2,3)ft 3 ' 4 > (2,1,3)} 

+25 3 (A; - k 2 ) [fi*> (1, 2, 3) - ^ (3, 1, 2)]^*(1, 2, 3) 

+25 3 (A; - k 3 )[^ 6 \2, 3, 1) - ft 6 \l, 2, 3)]^*(1, 3, 2)} . (A31) 




Figure 1: (Color online) The TMDs fx, g±L, h\ as functions of x and k\ are shown in the upper, 
middle and lower panels, respectively. Results for up and down quarks are given in the left (right) 
panels. 




Figure 2: (Color online) The TMDs H\l and h^ T as functions of x and k\ are shown in the upper 
and lower panels, respectively. Results for up and down quarks are given in the left (right) panels. 
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Figure 3: (Color online) The contribution to the TMD from the L z = ±1 wave components 
(left panel) and (right panel) from the L z = and L z = 2 wave components as a function of x and 
fcj_. Upper (lower) panels for up and down quarks. 
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Figure 4: The parton distribution hh${x) (left panels) and the transverse moment h^, ±9 (x) (right 
panels). Solid curves: results from the light-cone CQM model with the momentum wave function 



of Ref. 



42]. Dashed curves in the upper panels: results from the 



momentum wave function in the hypercentral model of Ref. 



panels: results from the spectator model of Ref. 



46, 



light-cone CQM model with the 



4j] ] . Dashed curves in the lower 



481 ] . Dotted curves: results from the bag model. 
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X X 



Figure 5: The transverse moments q (left panels) and g[j, q (right panels) as functions of x. 
Upper (lower) panels for up and down quarks. Solid curves refer to the result obtained with the 
light-cone CQM model; dashed curves refer to the Wandzura-Wilczek-type approximation. 
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Figure 6: (Color online) Transverse shape of the proton, /3relt(&±, n )//i(^j_)) assuming a struck 
u quark. The horizontal axis is the direction of Sj_ and n = S±, <j) n = 0. The shapes vary from 
the outer circle to the internal line as k± is increased from to 2.0 GeV in steps of 0.25 GeV. The 
left figure is the results when only the contribution from the P waves to h^ T is taken into account, 
and the right picture shows the total results when also the S and D waves are included. 




Figure 7: (Color online) The same as in Fig. [6] except that (fi n = ir. The shapes vary from the inner 
circle to the external line as k± is increased from to 2.0 GeV in steps of 0.25 GeV. 
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Figure 8: (Color online) Transverse shape of the proton, /3relt(&_l, n)/ f\{k\) assuming a struck 
d quark. The horizontal axis is the direction of S± and n = S±, 4> n = 0. The shapes vary from 
the inner circle to the external line as k± is increased from to 2.0 GeV in steps of 0.25 GeV. The 
left figure is the results when only the contribution from the P waves to h^ T is taken into account, 
and the right picture shows the total results when also the S and D waves are included. 




Figure 9: (Color online) the same as in Fig. [5] except that <j) n = tt. The shapes vary from the outer 
circle to the internal line as k± is increased from to 2.0 GeV in steps of 0.25 GeV. 
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